Search for: All records

Abstract In inland water covering lakes, reservoirs, and ponds, the gas exchange of slightly soluble gases such as carbon dioxide, dimethyl sulfide, methane, or oxygen across a clean and nearly flat air‐water interface is routinely described using a water‐side mean gas transfer velocity , where overline indicates time or ensemble averaging. The micro‐eddy surface renewal model predicts , where is the molecular Schmidt number, is the water kinematic viscosity, and is the waterside mean turbulent kinetic energy dissipation rate at or near the interface. While has been reported across a number of data sets, others report large scatter or variability around this value range. It is shown here that this scatter can be partly explained by high temporal variability in instantaneous around , a mechanism that was not previously considered. As the coefficient of variation in increases, must be adjusted by a multiplier that was derived from a log‐normal model for the probability density function of . Reported variations in with a macro‐scale Reynolds number can also be partly attributed to intermittency effects in . Such intermittency is characterized by the long‐range (i.e., power‐law decay) spatial auto‐correlation function of . That varies with a macro‐scale Reynolds number does not necessarily violate the micro‐eddy model. Instead, it points to a coordination between the macro‐ and micro‐scales arising from the transfer of energy across scales in the energy cascade. 

Abstract The Arctic–Boreal Zone is rapidly warming, impacting its large soil carbon stocks. Here we use a new compilation of terrestrial ecosystem CO₂fluxes, geospatial datasets and random forest models to show that although the Arctic–Boreal Zone was overall an increasing terrestrial CO₂sink from 2001 to 2020 (mean ± standard deviation in net ecosystem exchange, −548 ± 140 Tg C yr⁻¹; trend, −14 Tg C yr⁻¹;P < 0.001), more than 30% of the region was a net CO₂source. Tundra regions may have already started to function on average as CO₂sources, demonstrating a shift in carbon dynamics. When fire emissions are factored in, the increasing Arctic–Boreal Zone sink is no longer statistically significant (budget, −319 ± 140 Tg C yr⁻¹; trend, −9 Tg C yr⁻¹), and the permafrost region becomes CO₂neutral (budget, −24 ± 123 Tg C yr⁻¹; trend, −3 Tg C yr⁻¹), underscoring the importance of fire in this region. 

Abstract We examined the seasonality of photosynthesis in 46 evergreen needleleaf (evergreen needleleaf forests (ENF)) and deciduous broadleaf (deciduous broadleaf forests (DBF)) forests across North America and Eurasia. We quantified the onset and end (Start_GPPand End_GPP) of photosynthesis in spring and autumn based on the response of net ecosystem exchange of CO₂to sunlight. To test the hypothesis that snowmelt is required for photosynthesis to begin, these were compared with end of snowmelt derived from soil temperature. ENF forests achieved 10% of summer photosynthetic capacity ∼3 weeks before end of snowmelt, while DBF forests achieved that capacity ∼4 weeks afterward. DBF forests increased photosynthetic capacity in spring faster (1.95% d⁻¹) than ENF (1.10% d⁻¹), and their active season length (End_GPP–Start_GPP) was ∼50 days shorter. We hypothesized that warming has influenced timing of the photosynthesis season. We found minimal evidence for long‐term change in Start_GPP, End_GPP, or air temperature, but their interannual anomalies were significantly correlated. Warmer weather was associated with earlier Start_GPP(1.3–2.5 days °C⁻¹) or later End_GPP(1.5–1.8 days °C⁻¹, depending on forest type and month). Finally, we tested whether existing phenological models could predict Start_GPPand End_GPP. For ENF forests, air temperature‐ and daylength‐based models provided best predictions for Start_GPP, while a chilling‐degree‐day model was best for End_GPP. The root mean square errors (RMSE) between predicted and observed Start_GPPand End_GPPwere 11.7 and 11.3 days, respectively. For DBF forests, temperature‐ and daylength‐based models yielded the best results (RMSE 6.3 and 10.5 days). 

Abstract Tundra and boreal ecosystems encompass the northern circumpolar permafrost region and are experiencing rapid environmental change with important implications for the global carbon (C) budget. We analysed multi-decadal time series containing 302 annual estimates of carbon dioxide (CO₂) flux across 70 permafrost and non-permafrost ecosystems, and 672 estimates of summer CO₂flux across 181 ecosystems. We find an increase in the annual CO₂sink across non-permafrost ecosystems but not permafrost ecosystems, despite similar increases in summer uptake. Thus, recent non-growing-season CO₂losses have substantially impacted the CO₂balance of permafrost ecosystems. Furthermore, analysis of interannual variability reveals warmer summers amplify the C cycle (increase productivity and respiration) at putatively nitrogen-limited sites and at sites less reliant on summer precipitation for water use. Our findings suggest that water and nutrient availability will be important predictors of the C-cycle response of these ecosystems to future warming. 

Abstract Accounting for temporal changes in carbon dioxide (CO₂) effluxes from freshwaters remains a challenge for global and regional carbon budgets. Here, we synthesize 171 site-months of flux measurements of CO₂based on the eddy covariance method from 13 lakes and reservoirs in the Northern Hemisphere, and quantify dynamics at multiple temporal scales. We found pronounced sub-annual variability in CO₂flux at all sites. By accounting for diel variation, only 11% of site-months were net daily sinks of CO₂. Annual CO₂emissions had an average of 25% (range 3%–58%) interannual variation. Similar to studies on streams, nighttime emissions regularly exceeded daytime emissions. Biophysical regulations of CO₂flux variability were delineated through mutual information analysis. Sample analysis of CO₂fluxes indicate the importance of continuous measurements. Better characterization of short- and long-term variability is necessary to understand and improve detection of temporal changes of CO₂fluxes in response to natural and anthropogenic drivers. Our results indicate that existing global lake carbon budgets relying primarily on daytime measurements yield underestimates of net emissions. 

Abstract Wetlands are responsible for 20%–31% of global methane (CH₄) emissions and account for a large source of uncertainty in the global CH₄budget. Data‐driven upscaling of CH₄fluxes from eddy covariance measurements can provide new and independent bottom‐up estimates of wetland CH₄emissions. Here, we develop a six‐predictor random forest upscaling model (UpCH4), trained on 119 site‐years of eddy covariance CH₄flux data from 43 freshwater wetland sites in the FLUXNET‐CH4 Community Product. Network patterns in site‐level annual means and mean seasonal cycles of CH₄fluxes were reproduced accurately in tundra, boreal, and temperate regions (Nash‐Sutcliffe Efficiency ∼0.52–0.63 and 0.53). UpCH4 estimated annual global wetland CH₄emissions of 146 ± 43 TgCH₄ y⁻¹for 2001–2018 which agrees closely with current bottom‐up land surface models (102–181 TgCH₄ y⁻¹) and overlaps with top‐down atmospheric inversion models (155–200 TgCH₄ y⁻¹). However, UpCH4 diverged from both types of models in the spatial pattern and seasonal dynamics of tropical wetland emissions. We conclude that upscaling of eddy covariance CH₄fluxes has the potential to produce realistic extra‐tropical wetland CH₄emissions estimates which will improve with more flux data. To reduce uncertainty in upscaled estimates, researchers could prioritize new wetland flux sites along humid‐to‐arid tropical climate gradients, from major rainforest basins (Congo, Amazon, and SE Asia), into monsoon (Bangladesh and India) and savannah regions (African Sahel) and be paired with improved knowledge of wetland extent seasonal dynamics in these regions. The monthly wetland methane products gridded at 0.25° from UpCH4 are available via ORNL DAAC (https://doi.org/10.3334/ORNLDAAC/2253). 

Abstract In eddy covariance measurements, the storage flux represents the variation in time of the dry molar fraction of a given gas in the control volume representative of turbulent flux. Depending on the time scale considered, and on the height above ground of the measurements, it can either be a major component of the overall net ecosystem exchange or nearly negligible. Instrumental configuration and computational procedures must be optimized to measure this change at the time step used for the turbulent flux measurement. Three different configurations are suitable within the Integrated Carbon Observation System infrastructure for the storage flux determination: separate sampling, subsequent sampling and mixed sampling. These configurations have their own advantages and disadvantages, and must be carefully selected based on the specific features of the considered station. In this paper, guidelines about number and distribution of vertical and horizontal sampling points are given. Details about suitable instruments, sampling devices, and computational procedures for the quantification of the storage flux of different GHG gases are also provided. 

Abstract. Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO2 net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ETT), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr−1 for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr−1 for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×103 ± 0.9×103 km3 globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates. 

Warming-induced carbon loss through terrestrial ecosystem respiration ( Re ) is likely getting stronger in high latitudes and cold regions because of the more rapid warming and higher temperature sensitivity of Re ( Q 10 ). However, it is not known whether the spatial relationship between Q 10 and temperature also holds temporally under a future warmer climate. Here, we analyzed apparent Q 10 values derived from multiyear observations at 74 FLUXNET sites spanning diverse climates and biomes. We found warming-induced decline in Q 10 is stronger at colder regions than other locations, which is consistent with a meta-analysis of 54 field warming experiments across the globe. We predict future warming will shrink the global variability of Q 10 values to an average of 1.44 across the globe under a high emission trajectory (RCP 8.5) by the end of the century. Therefore, warming-induced carbon loss may be less than previously assumed because of Q 10 homogenization in a warming world.